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Abstract 

The state of a quantum system, consisting of two distinct subsystems, is called separable if 
it can be prepared by two distant experimenters who receive instructions from a common 
source, via classical communication channels. A necessary condition is derived and is 
shown to be more sensitive than Bell's inequality for detecting quantum inseparability. 
Moreover, collective tests of Bell's inequality (namely, tests that involve several composite 
systems simultaneously) may sometimes lead to a violation of Bell's inequality, even if 
the latter is satisfied when each composite system is tested separately. 
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1. Introduction 



Quantum entanglement is a physical resource. Applications include secure communication 

between distant observers [1], faithful teleportation of an unknown quantum state to an 
unknown location [2], quantum computation, and in particular error correction codes [3]. 
Quantum entanglement is subtler than classical correlation. A composite quantum system 
can be prepared in a prescribed state p, with distant correlations, by two observers who 
only use local raw materials and receive instructions from a common source by a classical 
communication channel. Such a preparation method yields a density matrix which is 
separable into a sum of direct products, 



where the positive weights wk satisfy J2wk = 1, and where and p'^ are density 
matrices for the two subsystems. A separable system represented by the above p is 
correlated, but it is not entangled: it always satisfies Bell's inequality [4]. 

That inequality was originally derived for solving a completely different problem: is 
quantum theory compatible with an underlying pseudo-classical "subquantum" back- 
ground (deterministic or possibly stochastic)? Such a post-quantum theory would pre- 
sumably involve additional "hidden" variables, and the statistical predictions of ordinary 
quantum theory would be reproduced by performing suitable averages over these hidden 
variables. Bell [4] was the first to show that if the constraint of locality is imposed on 
the hidden variables (namely, if the hidden variables of two distant quantum systems are 
themselves separable into two distinct subsets), then there is an upper bound to the cor- 
relations of results of measurements that can be performed on the two distant systems. 
That bound is violated by some states in quantum mechanics, for example the singlet 
state of two spin-| particles. It is not quantum theory itself which is nonlocal. It is 
any pseudo-classical theory which attempts to mimic quantum effects that is necessarily 
nonlocal and contextual. 

A variant of Bell's inequality, more general and more useful for experimental tests, was 
derived by Clauser, Horne, Shimony, and Holt (CHSH) [5] . It can be written 



On the left hand side, A and A' are two operators that can be measured by an observer, 
conventionally called Alice. These operators do not commute (so that Alice has to choose 
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whether to measure A or A') and each one is normahzed to unit norm (the norm of an 
operator is defined as the largest absolute value of any of its eigenvalues). Likewise, B and 
B' are two normalized noncommuting operators, any one of which can be measured by 
another, distant observer (Bob). Note that each one of the expectation values in Eq. (^ 
can be calculated by means of quantum theory, if the quantum state is known, and it is 
also experimentally observable, by repeating the measurements sufficiently many times, 
starting each time with identically prepared pairs of quantum systems. The validity of 
the CHSH inequality, for all combinations of measurements independently performed on 
both systems, is a necessary condition for the possible existence of a local hidden variable 
(LHV) model for the results of these measurements. It is not in general a sufficient 
condition, as will be shown below. 

Note that, in order to test Bell's inequality, the two distant observers independently 
measure subsytems of a composite quantum system, and then they report their results to 
a common site where that information is analyzed [6]. This is the converse of the situation 
that was mentioned above, where the two observers received instructions from a common 
center. There are density matrices for which it can be proved that no such set of in- 
structions exists, and yet Bell's inequality is satisfied [7-10]. I shall derive below a simple 
algebraic test which is a necessary condition for the existence of the decomposition (|I|). 
I shall then give some examples showing that this new criterion is more restrictive than 
Bell's inequality, or than the a-entropy inequality [11]. 



2. Separability of density matrices 

The derivation of the separability condition is easiest when the density matrix elements 
are written explicitly, with all their indices [6]. For example, Eq. (|ID becomes 

Pmfj.,nu — {Px)mn {,PK)tJ,u- (3) 

K 

Latin indices refer to the first subsystem, Greek indices to the second one (the sub- 
systems may have different dimensions). Note that this equation can always be satisfied 
if we replace the quantum density matrices by classical Liouville functions (and the dis- 
crete indices are replaced by canonical variables, p and q). The reason is that the only 
constraint that a Liouville function has to satisfy is being non-negative. On the other 
hand, we want quantum density matrices to have non- negative eigenvalues, rather than 
non-negative elements, and the latter condition is more difficult to satisfy. 
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It is noteworthy that any set of instructions that Ahce and Bob may receive is inherently 
ambiguous, if these observers have no way of ascertaining that what one of them calls 
i = v^— 1 is the same as i (not —i) of the other observer. Indeed, the only way of explaining 
to someone what i and —i actually are is by referring to a material (chiral) object, such 
as a screw. If we are restricted to the use of verbal instructions for preparing p, the result 
may be either p ot p* = . Therefore Alice and Bob may also end up with the "dual" 
matrix 

^=E^^(PAr®PA- (4) 

A 

This ambiguity is an important clue for solving the following problem. Let all the 
matrix elements of p be given. Is a decomposition such as in Eq. possible? A simple 
condition is readily obtained by defining a new matrix, 

C"m/i,ni' = Pnfj,,mv (5) 

The Latin indices of p have been transposed, but not the Greek ones. This is not a unitary 
transformation but, nevertheless, the a matrix is Hermitian. Whenever Eq. (|1]) is valid, 
the transposed matrices {p'^)'^ = {p'a)* non-negative matrices with unit trace, so that 
they also are legitimate density matrices. It follows that none of the eigenvalues of a is 
negative. This is a necessary condition for Eq. (|I|) to hold [12]. 

Note that the eigenvalues of a are invariant under separate unitary transformations, 
U' and U", of the bases used by the two observers. In such a case, p transforms as 

p^{U'®U")piU'®U'y, (6) 

and we then have 

a^{U'^(^U")aiU'^®Uy, (7) 

which also is a unitary transformation, leaving the eigenvalues of a invariant. 

As an example, consider a pair of spin-| particles in an impure singlet state, consisting 
of a singlet fraction x and a random fraction (1 — a;) [13]. Note that the "random fraction" 
(1— a;) also includes singlets, mixed in equal proportions with the three triplet components. 
We have 

Pm^,nu SjYifi^nv (l ^mn ^iiv /4, (8) 
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where the density matrix for a pure singlet is given by 

'S'01,01 = 'S'10,10 = — 'S'01,10 = — 'S'10,01 — |, (9) 

and all the other components of S vanish. (The indices and 1 refer to any two ortho- 
gonal states, such as "up" and "down.") A straightforward calculation shows that a has 
three eigenvalues equal to (1 -|-x)/4, and the fourth eigenvalue is (1 — 3x)/4. This lowest 
eigenvalue is positive \i x < |, and the separability criterion is then fulfilled. This result 
may be compared with other criteria: Bell's inequality holds for x < and the a- 

entropic inequality [11] for x < l/-\/3. These are therefore much weaker tests for detecting 
inseparability than the condition that was derived here. 

In this particular case, it happens that this necessary condition is also a sufficient one. 
It is indeed known that if x < | it is possible to write p as a mixture of unentangled 
product states [14]. This suggests that the necessary condition derived above {a has no 
negative eigenvalue) might also be sufficient for any p. A proof of this conjecture was 
indeed recently obtained [15] for composite systems having dimensions 2x2 and 2x3. 
However, for higher dimensions, the present necessary condition is not a sufficient one. 
Some counterexamples have been constructed, with dimensions 2x4 and 3x3 [16]. It was 
also shown [16] that the decomposition of any separable p requires at most (dim 7i)^ terms 
of rank one. This property is an important step toward finding an efficient algorithm for 
the decomposition of any separable p. Such an algorithm, which could also indicate that 
p is not separable, is still lacking at the time of writing. 

As a second example, consider a mixed state consisting of a fraction x of the pure state 
a|01) + 6|10) (with |a|^ + = 1), and fractions (1 — x)/2 of the pure states |00) and 



1 11). We have 

Poo,oo = Pii,ii = (1 - a;)/2, (10) 

Poi,oi ^ x\a\^ , (11) 

Pio,io = a^l&r, (12) 

Poi,io = P*w,oi = xab*, (13) 

and the other elements of p vanish. It is easily seen that the a matrix has a negative 
determinant, and thus a negative eigenvalue, when 

x> {l + 2\ab\)-\ (14) 



5 



This is a lower limit than the one for a violation of Bell's inequality, which requires [10] 

X > [I + 2\ab\{V2 - (15) 

An even more striking example is the mixture of a singlet and a maximally polarized 
pair: 

Pmii,nu = X Smij,,nu + (1 — x) 6mO ^nO ^fiO ^uO- (16) 

For any positive x, however small, this state is inseparable, because a has a negative 
eigenvalue (— x/2). On the other hand, the Horodecki criterion [17] gives a very generous 
domain to the validity of Bell's inequality: x < 0.8. 



3. Collective tests for nonlocality 

The weakness of Bell's inequality as a test for inseparability is due to the fact that the 
only use made of the density matrix p is for computing the probabilities of the various 
outcomes of tests that may be performed on the subsystems of a single composite system. 
On the other hand, an experimental verification of that inequality necessitates the use of 
many composite systems, all prepared in the same way. However, if many such systems 
are actually available, we may also test them collectively, for example two by two, or three 
by three, etc., rather than one by one. If we do that, we must use, instead of p (the density 
matrix of a single system), a new density matrix, which is p (S> p, or p p p, in a higher 
dimensional space. It will now be shown that there are some density matrices p that 
satisfy Bell's inequality, but for which p® p, oi p® p® p, etc., violate that inequality [18]. 

The example that will be discussed is that of the Werner states [7] defined by Eq. (|^). 
Let us consider n Werner pairs. Each one of the two observers has n particles (one 
from each pair). They proceed as follows. First, they subject their ra-particle systems to 
suitably chosen local unitary transformations, f/, for Alice, and V , for Bob. Then, they 
test whether each one of the particles labelled 2, 3, . . . , n, has spin up (for simplicity, it 
is assumed that all the particles are distinguishable, and can be labelled unambiguously). 
Note that any other test that they can perform is unitarily equivalent to the one for spins 
up, as this involves only a redefinition of the matrices U and V . If any one of the 2{n — 1) 
particles tested by Alice and Bob shows spin down, the experiment is considered to have 
failed, and the two observers must start again with n new Werner pairs. 
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A similar elimination of "bad" samples is also inherent to any experimental procedure 
where a failure of one of the detectors to fire is handled by discarding the results registered 
by all the other detectors: only when all the detectors fire are their results included in 
the statistics. This obviously requires an exchange of classical information between the 
observers. (There is a controversy on whether a violation of Bell's inequality with post- 
selected data [19] is a valid test for nonlocality [20]. I shall not discuss this issue here; I 
only examine whether or not Bell's inequality is violated by the postselected data.) 

The calculations shown below will refer to the case n = 3, for definiteness. The 
generalization to any other value of n is straightforward. Spinor indices, for a single spin-| 
particle, will take the values (for the "up" component of spin) and 1 (for the "down" 
component). The 16 components of the density matrix of a Werner pair, consisting of a 
singlet fraction x and a random fraction (1 — x), are, in the standard direct product basis: 

Pmn,st = xSmn,st + [l - x) 5ms ^nt / (17) 

where I am now using only Latin indices, contrary to what I did in Eq. (|^); this is because 
Greek indices will be needed for another purpose, as will be seen soon. Thus, now, the 
indices m and s refer to Alice's particle, and n and t to Bob's particle. 

When there are three Werner pairs, their combined density matrix is a direct product 
p®p'®p" , or explicitly, pmnM Pm'n',s't' Pm"n",s"t"- The result of the unitary transformations 
U and V is 

p ® p' ® p" ^ (f/ ® V) {p®p'® p") {U^ ® V^). (18) 
Explicitly, with all its indices, the U matrix satisfies the unitarity relation 

mm'm" 

In order to avoid any possible ambiguity, Greek indices (whose values are also and 1) 
are now used to label spinor components after the unitary transformations. Note that the 
indices without primes refer to the two particles of the first Werner pair (the only ones 
that are not tested for spin up) and the primed indices refer to all the other particles (that 
are tested for spin up). The K,,yV",nn'n" matrix elements of Bob's unitary transformation 
satisfy a relationship similar to (19). The generalization to a larger number of Werner 
pairs is obvious. 

After the execution of the unitary transformation (|l5), Alice and Bob have to test 
that all the particles, except those labelled by the first (unprimed) indices, have their 
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spin up. They discard any set of n Werner pairs where that test fails, even once. The 
density matrix for the remaining "successful" cases is thus obtained by retaining, on the 
right hand side of Eq. (p!8[), only the terms whose primed components are zeros, and then 
renormalizing the resulting matrix to unit trace. This means that only two of the 2" rows 
of the U matrix, namely those with indices 000. . . and 100. . . , are relevant (and likewise 
for the V matrix). The elimination of all the other rows greatly simplifies the problem of 
optimizing these matrices. We shall thus write, for brevity, 

where /i = 0, 1. Then, on the left hand side of Eq. (p!9[) , we effectively have two unknown 
row vectors, Uq and Ui, each one with 2" components (labelled by Latin indices mm'm"). 
These vectors have unit norm and are mutually orthogonal. Likewise, Bob has two vectors, 
Vo and Vi. The problem is to optimize these four vectors so as to make the expectation 
value of the Bell operator [21], 

C:=AB + AB' + A'B - A'B', (21) 
as large as possible. 

The optimization proceeds as follows. The new density matrix, for the pairs of spin-| 
particles that were not tested by Alice and Bob for spin up (that is, for the first pair in 
each set of n pairs), is 

(yPncw) fiu,aT = N U fj,^rnm' m" '^u,nn'n" Pmn,st Pm,'n',s't' Pm"n" ,s"t" U^^ggigii V^^tt't" y (^2) 

where iV is a normalization constant, needed to obtain unit trace (A^~^ is the probability 
that all the "spin up" tests were successful). We then have [17], for fixed pnew and all 
possible choices of C, 

max [Tr (Cpnow)] = 2y/M, (23) 

where M is the sum of the two largest eigenvalues of the real symmetric matrix T'^T, 
defined by 

Tpg ■= Tr [((Tp O O-g) Pnew]- (24) 

(In the last equation, cTp and are the Pauli spin matrices.) Our problem is to find the 
vectors 11^ and V^, that maximize M. 
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At this point, some simplifying assumptions are helpful. Since all matrix elements 
Pmn,st are real, it is reasonable to restrict the search to vectors and that have 
only real components. Furthermore, the situations seen by Alice and Bob are completely 
symmetric, except for the opposite signs in the standard expression for the singlet state: 

V'= [(;)(;)- fflOl (25) 



These signs can be made to become the same by redefining the basis, for example by 
representing the "down" state of Bob's particle by the symbol (^^^ , without changing the 
basis used for Alice's particle. This unilateral change of basis is equivalent a substitution 

Vu,nn'n" ^ { — ^Y'^"''^"' Vu,nn'n", (26) 

on Bob's side. The minus signs in Eq. (P) also disappear, and there is complete symmetry 
for the two observers. It is then plausible that, with the new basis, the optimal 1/^ and 
are the same. Therefore, when we return to the original basis and notations, they satisfy 

K/,nn'n" = ( — 1)^^"^" ^" Uu,nn'n" ■ (27) 

We shall henceforth restrict our search to pairs of vectors that satisfy this relation. 

After all the above simplifications, the problem that has to be solved is the following: 
find two mutually orthogonal unit vectors, Uq and Ui, each one with 2" real components, 
that maximize the value of M{U) defined by Eqs. (EHf) and (|2^. This is a standard opti- 



mization problem which can be solved numerically. Since the function M{U) is bounded, 
it has at least one maximum. It may, however, have more than one: there may be several 
distinct local maxima with different values. A numerical search leads to one of these 
maxima, but not necessarily to the largest one. The outcome may depend on the initial 
point of the search. It is therefore imperative to start from numerous randomly chosen 
points in order to ascertain, with reasonable confidence, that the largest maximum has 
indeed been found. 



4. Numerical results 

In all the cases that were examined, M[U) 
following simple choice: 

^0,00... = = 1, 



turned out to have a local maximum for the 

(28) 
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and all the other components of Uq and Ui vanish. Recall that the "vectors" Uq and Ui 
actually are two rows, f/ooo... and f/ioo..., of a unitary matrix of order 2" (the other rows 
are irrelevant because of the elimination of all the experiments in which a particle failed 
the spin-up test). In the case n = 2, one of the unitary matrices having the property (EBp 



is a simple permutation matrix that can be implemented by a "controUed-NOx" quantum 
gate [22] known as XOR (exclusive or). For larger values of n, matrices that satisfy 
Eq. ( 128|) will also be called xoR-transformations. 

It was found, by numerical calculations, that xoR-transformations always are the op- 
timal ones for n = 2. They are also optimal for n = 3 when the singlet fraction x is less 
than 0.57, and for n = 4 when x < 0.52. For larger values of x, more complicated forms 
of Uq and Ui give better results. For n = 3, it was recently found that the optimum is a 
"controlled Hadamard transform" 



^0,000 — f^o,iii — f^i.ooi — ^1,100 — (29) 

The corresponding result for higher n is not known. Anyway, the existence of two different 
sets of maxima may be seen in Fig. 1: there are discontinuities in the slopes of the graphs 
for n = 3 and 4, that occur at the values of x where the largest value of (C) jumps from 
one local maximum to another one. 

For n = 5, a complete determination of Uq and Ui requires the optimization of 64 pa- 
rameters subject to 3 constraints, more than my workstation could handle in a reasonable 
time. I therefore considered only xoR-transformations, which are likely to be optimal for 
a;<0.5. In particular, for x = 0.5 (the value that was used in Werner's original work [7]), 
the result is (C) = 2.0087, and the CHSH inequality is violated. This violation occurs in 
spite of the existence of an explicit LHV model that gives correct results if the Werner 
pairs are tested one by one. 

These results prompt a new question: can we get stronger inseparability criteria by 
considering p (g) p, or higher tensor products? It is easily seen that no further progress can 
be achieved in this way. If p is separable as in Eq. (|l|), so is p ® p. Moreover, the partly 
transposed matrix corresponding to p ® p simply is a ® a, so that if no eigenvalue of a is 
negative, then a ® a too has no negative eigenvalue. 
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Caption of figure 

FIG. 1. Maximal expectation value of the Bell operator, versus the singlet fraction in the 
Werner state, for collective tests performed on several Werner pairs (from bottom to top 
of the figure, 1, 2, 3, and 4 pairs, respectively). The CHSH inequality is violated when 
(C) > 2. 
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